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Liischer established a non-perturbative formula to extract the elastic scattering phases from two- 
particle energy spectrum in a torus using lattice simulations. Rummukainen and Gottlieb further 
extend it to the moving frame, which is devoted to the system of two identical particles. In this 
work, we generalize Rummukainen-Gottlieb's formula to the generic two-particle system where two 
particles are explicitly distinguishable, namely, the masses of the two particles are different. The 
finite size formula are achieved for both C± v and C2v symmetries. Our analytical results will be 
very helpful for the study of some resonances, such as kappa, vector kaon, and so on. 
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I. INTRODUCTION 

Many low energy hadrons, such as kappa, sigma, can 
observed as resonances in the experiments. The energy 
eigenvalues of two-particle systems can be achieved by 
calculating the propagators using lattice QCD. Hence, it 
is highly desired to connect these calculated energy eigen- 
states to the experimental scattering phases. Liischer p]- 
[j| have established a non-perturbative formula to connect 
two-particle system's energy in a torus with the clastic 
scattering phase. Rummukainen and Gottlieb further ex- 
tended Liischer's formula to the moving frame (MF) Q. 
Moreover, Xu Feng et al generalized Liischer's formula 
in an asymmetric box [3] . These formulae have been em- 
ployed in many different applications (cl-fl9| . 

For some cases, we have to use the generic two-particle 
system to extract the resonance parameters in the mov- 
ing frame. However, all of these aforementioned formu- 
lae in the moving frame can apply only to two identi- 
cal particle system. For example, to examine the be- 
havior of the k resonance, it is highly desired for us 
to investigate ttK scattering in the moving frame. For 
a generic two-particle system in the moving frame, we 
should change the Rummukainen-Gottlieb's formulae, 
which is devoted to the system of two identical particles. 
To this purpose, we will strictly derive the equivalents of 
the Rummukainen-Gottlieb's formulae for a generic two- 
particle system in the moving frame not only from the- 
oretical aspects, but also from practical considerations. 
This is very helpful for the lattice study since it provides 
a feasible method in the study of the n decay, vector kaon 
K* decay, and so on. 

The alterations of the Rummukainen-Gottlieb's formu- 
lae arc mainly relevant with the different symmetries of 
two-particle system in a torus . The representations 
of the rotational group 0{h) are decoupled into irrc- 
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ducible representations of the D^h and D^h cubic groups 
for the system of two identical particles with the non- 
zero total momentum in a torus @. In a generic two- 
particle system, the symmetry is reduced. For the case 
of d = (0, 0, 1), the basic group is C± v instead of D^; As 
for d = (0, 1, 1), the basic group becomes C% v . Hence, 
the final finite size formula for the generic two-particle 
system in the moving frame is certainly new. 

This paper is organized as follows. In Sec. |Hl we 
discuss the basic properties of the generic two-particle 
states in a torus. In Sec. HID we extend Rummukainen- 
Gottlieb's formulae to the generic case and derive the 
fundamental relationship for the scattering phase in 
Eq. and in Scc. lIVI we present the symmetry consid- 
erations. Finally we give our brief conclusions in Sec. [V] 



II. GENERIC TWO-PARTICLE SYSTEM ON A 
CUBIC BOX 

Here we derive the formalisms required for calculating 
the scattering phases in a cubic box, which are enough 
for studying lattice simulation data. However, in con- 
crete lattice calculation, we should address the lattice 
artifacts [2(| ■ In this section we follow the Rummukainen 
and Gottlieb's notations and conventions @ , and extend 
them to the generic two-particle systems. 

Without loss of generality, we consider two particles 
with masses mi and to 2 for particle 1 and particle 2, 
respectively. In this work we are specially interested in 
a system having a non-zero total momentum, namely, 
the moving frame @. Using a moving frame with total 
momentum P = (2w/L)d, d £ Z 3 , the energy eigenvalues 
for two-particle system in the non-interacting case are 
given by @ 



Emf — 



Pi 



where pi = |pi|, pa = IP2I, and pi, p 2 denote the three- 
momenta of the particle 1 and particle 2, respectively, 
which obey the periodic boundary condition, 
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and the relation 



Pi + P2 = P- 



(1) 



In the center-of-mass (CM) frame, the total CM mo- 
mentum disappears, namely, 



P 



P 2 > 



where p* = (2ir/L)n, and n £ 1? . Here and hereafter we 
define the CM momenta with an asterisk (*). The possi- 
ble energy eigenstates of the generic two-particle system 
are given by 



E, 



CM 



- P 



*2 



P 



*2 



The rclativistic four-momentum squared is invariant, 
and Ecm is related to Em f in the moving frame through 
the Lorcntz transformation 



TP 2 

^MF 



In the moving frame, the center-of-mass is moving with 
a velocity of v = P/Emf- Using the standard Lorentz 
transformation with a boost factor 7 = l/Vl — v 2 , the 
Ecm can be obtained through Ecm = 7~ 1 Pmf, and 
momenta and p* are related by the standard Lorentz 
transformation, 



Pi=7(p* + v£ x *; 



P2 



-7(p* - vE 



2) ' 



(2) 



where E^ and -EJ are energy eigenvalues of the particle 1 
and particle 2 in the center-of-mass frame, respectively, 



Ei 



1 



2E CM 
1 (E? 



Ecm 



2E C 



M 



CM 



(3) 



and the boost factor acts in the direction of v, here and 
hereafter we adopt the shorthand notation 



7P = 7P|| + p ± , 



7 x p = 7 



(4) 



where p|| and pj_ are the p components which are par- 
allel and perpendicular to the CM velocity, respectively, 
namely, 



Pll = 



p v 



P - Pll 



(5) 



Thus, by inspecting Eqs. ([T]), and ©, it can be seen 
that the p* are quantized to the values 



where the set Pd is 



Pa 



2tt 
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r€Pd, 



E 2 

^CM 



(G) 



In the interacting case, the Ecm is given by 



Ecm = + k 2 



+ k\ k = ^q. (7) 



where q is no longer required to be an integer. Solving 
this equation for scattering momentum k we get 

1 r^o. ~2 

k = -= v [Ecm - ( m i- m2 ) 2 ]i E CM - Oi + m 2 ) 2 ]. 

£&CM 

(8) 

We can rewrite Eq. (|8]) to an elegant form as 

= 2 



E 



CM 



4? 



£E 2 



(9) 



CM 



It is exactly this energy shift between the non- 
interacting situation and interacting case, namely, 
Ecm — Ecm (or cquivalently |n| 2 — q 2 ), that we can 
evaluate two particle scattering phase. 

As it is done in Ref. @ , in the current study, we mainly 
investigate two moving frames. One is d = (0,0,1), 
where energy eigenstates transform under the tetragonal 
group Ca v , only the irreducible representation A\ is rele- 
vant for two-particle scattering states in a cubic box with 
angular momentum / = 0. Another one is d — (0, 1, 1), 
where energy eigenstates transform under the tetragonal 
group C2V, only the irreducible representation A\ is rel- 
evant. For the other cases, like d = (1, 1, 1), etc., we can 
easily work out from same way without difficulty. 

Assuming that the phase shifts Si with I = 1,2,3,... 
are tiny in our concerned energy range, the s-wave phase 
shift 5q is connected to the scattering momentum k by 



77T 3 / 2 <J 



tanJo(fe) = -^dTT 



(10) 



where k = (2tt /L)q 1 and the modified zeta function is 
defined by 

1 



r£P d 



(\r\ 2 -q 2 y 



(11) 



and the set Pd is 
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r = 7 1 


rH 


-H 



p2 

^CM 



n G 



(12) 

For Eq. (fTU|) . we note that the almost same result has 
already existed in Eq. (1) of Ref. [2l|, where the for- 
mula was just presented without any explanation. We 
can view our work as further confirming and strictly 
proving this formula. The modified zeta function con- 
verges when Re 2s > I + 3, and it can be analytically 
continued to whole complex plane. The scattering mo- 
mentum k is defined from the invariant mass ^/s as 
y/s = Ecm = \fk 2 + m\ + y k' 2 + m\ . The calculation 
method of Z^ {l]q 2 ) is discussed in Appendix IA1 and in 
Ref. [2(j. Using Eq. (fTU|) . we can obtain the phase shift 
from the energy spectrum using lattice simulations. If we 
now set ttl\ — m 2 , all the results in Ref. Q are elegantly 
recovered. 
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III. DERIVATION OF THE PHASE SHIFT 
FORMULA 

Here we deduce the basic phase shift formula in 
Eq. (|10|) for the generic two-particle system with spin- 
0. We utilize the Rummukainen and Gottlieb's formulae 
in Ref. @, and extend them to the generic two-particle 
system. To make the derivation simple, we are studying 
the system by the rclativistic quantum mechanics. 

Throughout this section, we employ the metric ten- 
sor sign convention = diag(l, —1, —1, —1), write the 
scalar productions in a compact form p 2 = p ■ p = p^p^ , 
etc, and express the quantities in natural units with 
H = c = 1. Here and hereafter we follow the original 
notations in Refs. ral. 



where 



P 



m 2 pi - mip 2 



mi + m 2 

P = pi + p 2 , 
M = Till + J77.2, 



(19) 

(20) 
(21) 



p is relative four-momentum operator, P is total four- 
momentum operator, and M is total mass of two parti- 
cles. Adding 1/mi x (fTT)) to l/m^x (fT5)) and subtracting 
(fT7|) from (fl"8j) . respectively, yield 



M 2 



mi — m 2 

2M 



M 2 



P^ 



i>(x,X) = 0, (22) 
tp(x,X) = 0. (23) 



A. Lorentz transformation of wave function 

Let us consider the generic system of two spin-0 parti- 
cles with mass m 1; and m 2 , respectively, in a cubic box. 
The two-particle system is described by the scalar wave 
function ip(xi,X2), where Xi = (a;^,x,), % = 1,2 are the 
four-dimensional Minkowski coordinates of two particles. 
The wave function changes under the standard Lorentz 
transformations, namely, 



ip(xi,x 2 ) = ip' {x'i,x 2 ) = ip' (Axi, Ax 2 ), 



(13) 



where (x'Y = K^ v x v defines the standard Lorentz trans- 
formation of the four-vector x. 

We can make the problem easier through the special 
properties of the CM frame. First let us study two non- 
interacting particles, and the wave functions of the sys- 
tem obey the Klein-Gordon (KG) equations 



(pi^Pi 11 ~ rnf) ip(xi,x 2 ) = 0, 
(hfj.pY - ml) i/j(x 1 ,x 2 ) = 0, 



(14) 



where pi^, i = 1,2 is the four- momentum operator. It 
is well-known, the problem simplifies if we separate the 
variables under the transformations 



X = 



mix\ + m 2 x 2 

m\ + m 2 
xi - x 2 , 



(15) 
(16) 



where X is the position of the CM, and x is the relative 
coordinate of two particles. Let us restrict ourselves to 
the solutions which are the eigenstates of the CM mo- 
mentum operator. Then Eq. (|14[) can be changed into 



Jf 2 P,P"+P,P 



.p» + m 2 



^P^+p^ + ^-P^ + m 2 



i>(x,X) =0,(17) 
1>{x,X) =0,(18) 



It is well-known that, without external potentials, the 
total momentum of the two-particle system is conserved, 
thus we can restraint ourselves to the eigenfunctions of 
P, namely, 



i/>(x, X) = e 



<j>{x), 



(24) 



where P M is a constant time-like vector, and P is denoted 
through P 2 =P fl P> 1 . 

In the present study, we are specially interested in the 
CM frame, which is denoted as the frame without the 
spatial components of the total momentum for the sys- 
tem, namely, P* = 0. Thus, we can only take the positive 
kinetic energy solutions P * = Eq M > mi + m 2 into con- 
sideration. So, Eqs. (|2"2"j) and (|2U)) can be rewritten as, 



PuP' 



^CM m l m 2 

(mi + m 2 ) 2 



Pa 



Ecm mi - m 2 
2 mi + m 2 



mim 2 4>cm(x*) 



<I>Cm(x*) 



275cm 



0, (25) 
0. (26) 



Eq. (f2l))) indicates Pq4>cm(x*) ^ for mi ^ m 2 . By in- 
specting Eq. (|23|) and Eq. ([2"B]) , we can reasonably assume 
that the wave function 0cm (x*) can be expressed as 



(j>CM(x*) = e lf3x 0cm(x*), 



(27) 



where x*° = xl° — x 2 ° is the relative temporal separation 
of two particles, and j3 is a constant, namely, 



Ecm m 2 - mi 
2 mi + m 2 



mj 



m\ 



2 Eqm 



(28) 



It is obvious that when mi = m 2 , (3 — >• 0. So, in the CM 
frame the wave function depends explicitly on the time 
variable t* = X*° 
relative spatial separation x* = x^ 
temporal separation of the particles x* , namely 



miXi° + m 2 x 2 °) I (mi + 7712)5 the 
xj, and the relative 



i>cu(x*,t*) = e - !BcMt *e^*°0 CM (x*), 



where the constant j3 is denoted in Eq. 



(29) 
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Now let us discuss the case in the moving frame. The 
transformation from the MF to CM frame can be ex- 
pressed as r* M = A^ u r v , where r is any position four- 
vector and quantities without * stand for these in the 
moving frame. Using the notation in (j4]), we have 



7(7" + v • r), r* = 7 (r + vr ), 



(30) 



where 7 is a boost factor, and v = P/Po is the three- 
velocity of the CM in the moving frame. We can rewrite 
v to a form for later use as 



2n 



2tt 



LEmf jLEcm 



(31) 



Considering the identity P^X^ = P*X*^, the Lorentz 
transformation in Eq. (fT3")) , and Eq. (|24"|) . the wave func- 
tion in the moving frame can be expressed as 



This result is consistent with the solution in Rcf. [22j . 

The Eqs. (f3"4"|) and (f3"T|) are very important when we 
study the wave functions of our system. Thus Eq. (|37l) is 
a very important result, which represents one of the main 
results of the present work. In the following studies, we 
take away the superscript * from the quantities in the CM 
frame. We can easily check that if we take toi = rri2 , all 
the corresponding results in Ref. [f| are restored. 



B. Modified singular d-periodic solutions of the 
Helmholtz equation 

In our concrete problem, for the potential V^(x) with 
a limit range @ , namely, 



U M (x)=0 for|x|>i?, 



(39) 



4>mf(x), 



where 



(/)mf(x) = ())mf(x , x) = (f> C M (7(2=° + v • x), 7 (x + vx )) . 

(32) 

Thus, the wave function </>mf depends on time separa- 
tion x° — x\ — x® explicitly. However, in the moving 
frame we only consider two particles with the same time 
coordinate, namely, a; = 0. It corresponds to the tilted 
plane (x*°, x*) = (7V • x, 7 x) for the CM frame, since the 
wave function </>cm is dependent of the relative temporal 
separation x*°, we can clearly observe the effect of the 
tilt to the wave function, and Eq. (pj2"j) take as 



4>MF (0, x) = (f) C M (7 V • x, 7 x) . 



(33) 



Using Eq. (|27|) and Eq. (j3Tj), we can rewrite Eq. (J33J) as 

0mf(O, x) = e ^ d - x / L CM (7x). (34) 
where j3' is a constant, namely, 



ra-2 — mi m\ — m\ 



mi + m 2 



p2 



(35) 



Eq. (|34p has a simple physical meaning: the CM system 
watches the torus in the moving frame expanded by a 
boost factor 7 in the direction of P, and the length scales 
in perpendicular directions are hold. At last, Eq. ([34[) 
connects the MF wave function. 



?/W(0,x,t,X) = e 



-iE MF t+iP-X 



</>mf(0,x), 



(36) 



to the CM frame wave function Eq. (|29p . By inspecting 
Eqs. ([25[) . (|26|) and ([29)) , we, at last, achieve the wave 
function (j)CM meeting the Helmholtz equation (HE) 



k ) 



'K'U 



(x*)=0, 



(37) 



where 



we suppose that the KG equation ([14]) in the CM frame 
still has a square integral solution. In the CM frame the 
interaction of the system is spherically symmetric. The 
wave function of the system is usually given in spherical 
harmonics 



0CM (x) 



E 



Yi m {B,(p)4>i m {x) 



(40) 



For x > R, </>cm is a solution of HE, and the radial 
functions </>; m meet the differential equation 



da; 2 



2_d_ 

x dx 



1(1 + 1) 



k 2 



4>lm(x) = 0, 



where 



k 2 = 



p2 



[m 



4E CM 



(41) 



(42) 



With the linear combinations of the spherical Bessel func- 
tions the solutions of Eq. (|4"Tj) can be expressed as 



4>im(x) = ci m [ai(k)ji(kx) + bi(k)ni{kx)} 



(43) 



Although we do not know the radial equation in the 
region x < R, by comparing the wave functions denoted 
in Eqs. ([40I43[) , we can get the relation between the scat- 
tering phase and the coefficients ai and bi Q, namely, 

e i2S t (k) = ai(k) + ibj(k) 
ai(k) - ibi[k) 

Because the radial equation ai and bi can be taken to be 
the real number for k > 0, Si(k) is real. Thus, for a given 
Z-sector, the phase shift 5i(k) can be expressed by energy 
-Emf through 



fi MF 



2\2 



1 (mf — m 2 ) 
4 E MF 



m\ 



ml 



1*2 



E 2 CM , (m 2 - ml) 2 mj + to. 



4^m 



(38) 



In the moving frame, we now investigate two generic 
particles confined in a torus of size I? with periodic 



5 



boundary conditions(PBC). The temporal direction of 
the torus is chosen to be infinite. The wave functions 
tpMF in MF should be periodic with regard to the posi- 
tion of each particle, namely, 

V , mf(xi,x 2 ) = V>mf(xi + li,X2 + mi), 1, m G I?. 

(44) 

The wave function i/>mf is provided by (J36J) , namely 



V>mf(xi,x 2 ) = exp 



P • (mixi + m 2 x 2 ) 

mi + 7712 



Combining Eq. (|4"4"|) and Eq. (|4"5")) gives 
2tt . 



W(xi-x 2 ). 
(45) 

(46) 



(x) = e 



'MF 



(x + nL), (47) 



where n = 1 m, d, n € Z 3 and P is the total mo- 
mentum 2 . The rule (|47| separates the wave functions 
into the different sectors, which can be categorized by 
the three- vector d. In the current study, we naturally 
concentrate on sectors d = (0, 0, 1) and d = (0, 1, 1). 

We are now on the position to employ Eq. (fMf to get 
the corresponding periodic rule for the CM wave func- 
tion. For a chosen vector d, we have 



^cm(x) = e 



i7rd-n ( 1 + 



6cm(x + 71iL), neZ 3 . 

(48) 

For simplcncss, we refer to the functions obeying the pe- 
riodicity rule (|4"8j) as modified d-periodic rule. As shown 
later, the modified d-periodic rule (|48]l is a milestone in 
this work. 

Now we deduce the general solutions of the HE satisfy- 
ing the modified d-pcriodic rule (|4"8"|) . Except the modi- 
fied d-pcriodicity, our derivation follows the works in sec- 
tion 4.4 of Ref. [f| . In this work, we will refer to a function 
as a singular modified d-periodic solution of the HE, if 
it is a smooth function denoted for all x ^ 7nL, n € Z 3 , 
and it meets the HE, namely, 

(V 2 + fc 2 )^(x) = 0, 

for k > 0, and meets the modified d-periodic rule, i.e., 

0(x) = e imrd "V(x + 7nL), n e Z 3 , (49) 

here and hereafter, for compactness, we defined a factor 
a, namely, 



a = 1 



p2 



(50) 



Eq. l|47|l can also be written as 



</>mf(x) 



»i+m 2 d ' n MP ( x + nL) . 



Following the almost same procedures and addressing the cor- 
responding formulae, we can arrive at the same final numerical 
results. For our case, d • n is an integer. 



When mi = m2, a — 1, this is the case Rummukainen 
and Gottlieb studied in Ref. @. Moreover, the wave 
function is required to be bounded by a power of l/|x| 
around the origin, namely 



lim |x A+1 <; 



< oo 



for the positive integer A, the degree of 4>. For our aim, 
it is enough to study the regular values of k, namely 



k + 



2- 

T 



7 



a 



We denote the Green function 
G d (x;fc)= 7 



per 



n G Z 3 



»tpx 



p 2 — k 2 



(51) 



where 

r = 



p e 



2tt 



n G T- 



Eq. ([51]) is well-defined due to the non-singular k. If we 
select k = (2n/L)^f- 1 (m + f d) , m G Z 3 , then 

k ■ (x + jnL) = k ■ x + and ■ n + 27rm • n, 

where n G Z 3 , and the Green function G d (x; k) meets 
clearly the modified d-pcriodic rule, as we expected. Fur- 
thermore, it obeys 



(V 2 + fc 2 )G d (x; k) = - 



e 

nGZ 3 



5(x + 7 nL) 



We can easily verify that the function G d (x; k) is a sin- 
gular modified d-periodic solution of Hclmholtz equation 
with degree 1. Other singular periodic solutions can be 
easily achieved by differentiating G d with x. And 



G d m (x;fc)=;y ;m (V)G d (x;fc), 



(52) 



where the harmonic polynomials iV; m (x) = x l Yi m (9,ip). 
We can easily verify that the functions Gf m arc singu- 
lar modified d-periodic solutions of the HE, and they 
make up a complete set of solutions so that any singular 
modified d-periodic solution of degree A is just a linear 
combination of Gf m (x;p) for I < A Q. In the region 
< x < L/2 the functions Gf m can be expanded in usual 
spherical harmonics, namely 



(-l)'fc' +1 



47T 

oo / 

E E 



[m(kx)Yi m (6, <p) 



ML l , m ,(k)j v (kx)Y l , m ,(e,<p) 



.(53) 



The regular part has the coefficients Aif m Vm ,{k). In the 
concrete calculation, only a few of them is needed, and 
we also provide its general expression, namely, 



,{k) 



i-iy 



3/2 



E E^t^ 2 n 



j=\l-l'\s= 



m,js,V m' i 



(54) 
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where q = kL/ '(2tt). Using Wigner 3j-symbols [23| the 
tensor Cj m j s< i' m > can be expressed as 



m.js.l'm' 



(_!)"» v /(2/ + l)(2j + l)(2i' + l) 

i j z' \ / Z j Z' 



m s —m 





The modified zeta function is formally denoted by 

yirn(r) 



r£P d 



(r 2 -g 2 ) s ' 



(55) 



(56) 



where 
Pd 



r e R 3 



7 



+ |d), neZ 3 } 



In Table U we summarized the expressions of Aif m v 
for I, V < 3. For compactness, we denoted 



»7, 



7r 3 / 2 v / 2T+T 



"^-'-^^(ljg 2 ). (57) 



The Wigner 3j-symbol values can be obtained in 
Ref. [iH ■ We can easily verify that, if we set mi = m 2 , 
all of the above definitions and formulae nicely reduce to 
the those obtained in Ref. @, as we expected. Of course, 
if we further select d = 0, the moving frame and the CM 
frame coincide, 7 — > 1 and Pd — > Z 3 , and they further 
neatly reduce to that in Ref. [|[ . The Table Q] can be 
compared with the Table 3 in Ref. @ , which summaries 
the matrix elements for mi = m 2 . We should stress that 
the functions wio, W30 and W50 appear in Tabic [I] If we 
set mi = TO2, then w\q — > 0, W30 — > 0, and w^o — > 0, and 
Rummukaincn-Gottlieb's results is elegantly restored. 
So far. we can arrive at 



E E vi m Gf m (x,k 2 ) = 

1=0 m=-l 
A I 

E E cim[ai(k)ji(kx)+bt(k)ni(kx)]Yi m (0,(p) (58) 



1=0 m=-l 

for the constants c; m and vi m . And we obtain 



A 



i' 



Q m ai( fc ) = E E C i'™' 6 ''( fc )- A/ Jrm',im( fc )- (59) 
l'=0 m'=-V 

■M-i'm' ,im can be viewed as the matrix element of the op- 
erator M. We can express Eq. (|5T))) as a matrix equa- 
tion, C(A — BM) = 0, where matrix An m \ni m i) = 
a i(p)5i t i>5 m ,m/ (similar for B). We can denote the phase 
shift matrix [D, 0] , 



e 2ii = 



A + iB 



A - iB 

The determinant condition arrives at 



det [e 2lS (M - i) - (M + i)] = 0. 



(60) 



IV. SYMMETRY DISCUSSIONS 

When the moving and CM frames coincide, the two- 
particle system has a cubic symmetry and the wave func- 
tions transforms under the representations of the cubic 
group Oh- On the other hand, according to Eq. ([34]) . 
if two frames are not equivalent, the Lorentz translation 
boost from the moving frame to the CM frame and only 
some subgroups of Oh group survive 0. 

In this work, we are mainly interested in a boost along 
one of the coordinate axes, namely d = (0, 0, 1). The ge- 
ometry of the torus transforms as (1,1,1)— > (1,1,7), an d 
the corresponding symmetry group is tetragonal point 
group Ci v , which has 8 elements: 4 rotations through 
an angle (wr/2), where n = 0, 1, 2, 3, around the £3-axis; 
and all four of the above multiplied by the reflection with 
respect to the (1, 3) -plane. The relevant point groups and 
the boost vectors are classified in Table HU 

In this work we mainly discuss the cubic and tetrag- 
onal symmetry groups Oh, C^ v , and C^. The tetrago- 
nal group Civ has four 1-dimcnsional representations A\ , 
Ai, B\, Bi, and one 2-dimensional representation E [24| . 
The representations of the rotational group are reduced 
into irreducible representations of C\ v as: 



r(°) = a x , 
r (1 > = Ai®e, 

T<- 2) = A 1 ®B 1 ®B 2 ®E. 



(61) 



The representations can be obtained through the char- 
acter tables [2^ or by enumerating harmonic polynomials 
of degree I which transform under the representations of 
Ci v ■ The basis polynomials for the corresponding repre- 
sentations are summarized in Table HTT1 for / < 2. 

The tetragonal group C% v has four 1-dimcnsional rep- 
resentations A\ , A2, B\, B2 [2H . The representations of 
the rotational group are further reduced into irreducible 
representations of C 2v as: 



r<°> = A lt 

r (1) = A 1 ®B 1 ®B 2 , 
r (2) = A 1 ®A 2 ®B 1 ~ 



(62) 



Bo 



Similarly, we can obtain the basis polynomials for C 2v 
representation, which are listed in Tabic ITVl for / < 2. 

In a typical lattice calculation, we can easily investigate 
is the Ai sector. We will therefore concentrate on this 
sector. As is seen, up to I < 2, s-wave, p-wave and d-wave 
contribute to this sector. The other symmetry sectors can 
be easily worked out in the same way. 

First, let us consider the angular momentum cutoff 
A = 0. From the reduction relations (|61I62[) and Ta- 
bles IIII1 IIV1 we note that only A4q Q 00 belongs to this 
sector, and Eq. (|6U)) can be written as 



tan 5o(k) 



1 



7<77T 



3/2 



- /Vl 00,00 



2 d (i;g 2 ) : 



2jt 



k. (63) 



This is our essential result. 
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I m I' m! 



lvl lm,l'm' 







"'00 



1 
1 

1 1 





1 

1 1 



woo 
woo 



+2w 2 o 

-W 2 



Wio 
WlO 



WOO 

woo 



woo 



+ f W20 




+ y-W40 
-y-W4.0 
|\/70W44 
+ f TO40 



-11030 



-jVUw 2 Q 



— jV21W40 
+ §\/l4u>40 



2\/3w44 



-Vis 



WlO 




-TV 


ifjWsO 


-*fv 


TTT W 50 




Jtt W5 ° 



WOO 
WOO 



"'00 



woo 



\--gW20 
+W20 



'3^20 



+ TJW W 
+ ^-W 4 



70W44 



7W40 



+ -33-W60 
-ffweo 



^-V42to 4 4 



14w 6 4 



-^V210to 6 4 



+— W40 



TABLE I: The matrix elements Mf ml , m , for d = (0, 0, 1), mi m 2 and I, I' < 3. 



d 


point group 


classification 


N elements 


(0,0,0) 


Oh 


cubic 


48 


(0,0, a) 


C4U 


tetragonal 


8 


(0,a,a) 


C2„ 


orthorhombic 


4 



TABLE II: The categorizations of the boosts on a cubic box 
and its reduction of the cubic group. The first column shows 
the boost direction and a is a non-zero real number. Here we 
adopted the Schonflies notation [24|. 



representation 


I = 


I = 1 


I = 2 


indices 


Ai 


1 


X3 


2 12 




A 2 










Bi 






T 2 T 2 

X\ — x 2 




B 2 






X1X2 




E 




Xi 


XiX-i 


i = 1,2 



TABLE III: The basis polynomials for the irreducible repre- 
sentations of Ci v . 



If the angular momentum cutoff A = 1, we should 
take the sector I = 1 into consideration, and the matrix 
in Eq. (|60p is two-dimensional. Hence, the determinant 
condition contains both phase shifts So and 5\, 



,2iS 



(m o 
10 (' 



m 



-i)- 

,2i<5 



(moo+i)] [e 2Ml (w*n 



1 



t> 2t<5i 



1 



i) - {mn + 0] 
(64) 



where we denote m a h = -M^ 60 . If 5i = (mod 7r), as 
what we expected, Eq. (|64|) reduces nicely to Eq. (|63| . 
If 5i 7^ 0. Normally we can reasonably suppose that the 
lowest i-channcl dominate the scattering phase. This is 



quite right in the low-energy elastic scattering 0, Q . 

If we expand So = Sq + Aq, where 8q meets Eq. (f6"5)l 
and Ao is a perturbative term, we can give the first order 
correction, namely, 

A (k) = a(k)5i(k). 

The function a(k) stands for the sensitivity of the higher 
scattering phases. For C& v symmetry, it is given by 



a(k) = - 



no 



m lo 



1' 



8 



representation 


I = 


I = 1 


I = 2 


Ai 


1 


X3 


x 3 ~ 


A 2 






a;ia;2 


Si 




Xl 


a;ia;3 


B 2 




X2 


2:2^3 



TABLE IV: The basis polynomials for the irreducible repre- 
sentations of Ci v . 



which is not usually small. To make sure that the 
Eq. (|63|) is a well approximation, the phase shift 
(5i(fc)(mod7r) should be small. Fortunately, the physical 
case is often like this way: the elastic scattering is often 
dominated by the lowest I- channel. 

The sensitivity function u(q 2 ) can be calculated using 
the matrix elements provided in Eq. (|57[) discussed in Ap- 
pendices and [B] The sensitivity function cr(q 2 ) for C± v 
symmetry with a = 1.15 and 7 = 1.177 is illustrated in 
Fig. [U here 7 is a boost factor, and a-factor is defined in 
Eq. (fST))) . which are typical values we used in Ref. [20(. In 
Fig. [U the lower panel is the same sensitivity function as 
the upper panel just to display some detailed variations. 
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v . , \J 




2 3 
q 2 



that it often varies in the range — 20, in Figs. [2] |3l 
and 21 we plotted just three of them. We can note that 
the sensitivity function <j{q 2 ) is finite for all q 2 > 0. For 
some special values of q 2 , however, the sensitivity func- 
tion a(q 2 ) has a sharp peak. For other values of q 2 away 
from these values, <j{p 2 ) is always not too large. 

We also note that, when q — > 0, the sensitivity function 
cr(g 2 ) is often large. However, it does not usually cause 
any problem because it is nicely canceled out by <5i which 
is of q 3 order at small q 0, Q. Thus, for the range < 
q 2 < 1.1 (except some special q 2 values), the Eq. (|63|) 
can be considered a good approximation. In fact, it is 
the range which are usually used to study the elastic 
scattering from lattice QCD [20]. 
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FIG. 2: The sensitivity <r(q 2 ) for C^ v symmetry with param- 
eters a = 1.05 and 7 = 1.177. 
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FIG. 3: The sensitivity a(q ) for C± v symmetry with param- 
eters a = 1.1 and 7 = 1.067. 



FIG. 1: The sensitivity cr(q 2 ) for Ca v symmetry with param- 
eters a = 1.15 and 7 = 1.177. 

In the work, we also calculate the sensitivity function 
o{q 2 ) using the some typical a and 7 values, we found 



We should bear in mind that if 8\{k) is not small, it 
is very difficult to extract the phase shift from energy 
spectrum Q. In principle, we can still extract the s- 
wave scattering phase from Eq. (|64[) through dividing the 
p-wave phase shift by lattice simulations at various en- 
ergy Q , because the corrections due to higher scattering 
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standing the dynamical aspect of hadron reactions with 
QCD. Moreover, BES collaboration recently carried out 
some experimental measurements pBl I2H ] to investigate 
k resonance mass and its decay width. With the modi- 
fied formula in Eq. (|63[) and our strict discussion of this 
formula from theoretical aspects, now it will be possible 
to compute the resonance masses and perhaps its de- 
cay widths of some resonances including possible exotic 
hadrons as well as traditional hadrons like k and vector 
kaon K* , etc., directly from lattice simulation in a correct 
manner. We have already used these formulae to prelimi- 
narily analyze our irK scattering at / = 1/2 channel [20j |. 
and the reasonable results of our lattice simulation data 
supports these formula. 



FIG. 4: The sensitivity cr{q 2 ) for Cav symmetry with param- 
eters a = 1.15 and 7 = 1.067. 
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phases can also be estimated from lattice calculations. 
For example, from Table IIII1 it is easily noted that, for 
lattices with C^ v symmetry, by inspecting energy eigen- 
values of E symmetry, we can get an approximate esti- 
mate for the p-wave scattering phase Si which dominates 
this symmetry sector. It seems to be too difficult, how- 
ever naturally, it is still possible to calculate the energy 
spectrum, this is our future tasks. 

If we choose the sector d = 0, the moving and 
the CM frames coincide, 7 — > 1 and — > Z 3 , 
(|64p nicely reduces to the form given in 
Of course, if we select mi = m,2 and Pd — > 
|r = 7- 1 (n + d/2)}, n g Z 3 , and Eq. l|B3|l 
neatly reduces to the form presented in Ref. @. These 
are what we expected. 

As for A = 2 or higher, it is quite complicated. See 
the relevant discussions in Ref. Q ■ Bearing in mind that 
this work is an exploratory study for some systems like 
the 7rA" system, the main purpose is to present some con- 
ceptual and theoretical issues. 
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Appendix A: The calculation of zeta function 

The method for evaluating the zeta function for d = 
has been discussed by Luscher in Ref. 0] . Rummukainen 
and Gottlieb extended this discussion in the MF for 
d^0,Q = 1 @. The formalism used here is further 
adapted to the case of d 7^ 0, a 7^ 1, and we just present 
the essential formulae. 

We first denotes the heat kernel on a modified d- 
periodic torus in Eq. (|4"5|). namely, 



— Y 



(Al) 



where the summation for r is carried out over the set 

a 



Pd = {1 



7 



-1 



n 



,;d), neZ 3 } 



(A2) 



In the current work we have strictly investigated 
the scattering states of two-particle with unequal mass, 
and the best-efforts are paid to derive the modified d- 
pcriodic rule which is crucial to the alteration of the 
Rummukaincn-Gottlicb's formula. The finite size expres- 
sions, which can be regarded as a generalization of the 
Rummukaincn-Gottlicb's formulae to the generic two- 
particle system in the moving frame, are developed. We 
also checked that all the Rummukainen-Gottlieb's results 
in Ref. 0] is nicely restored if we set mi = mi. 

Since the so-called k meson is a low-lying scalar meson 
with strangeness, a study of k meson decay is an explicit 
exploration of the three-flavor structure of the low-energy 
hadronic interactions, which is not directly probed in mr 
scattering, therefore, it is a significant step for us undcr- 



here the factor a is denoted in Eq. ([50]) . and the operation 
7 _1 is is defined in Eq. @. Following from Poisson's 
identity, we can rewrite the heat kernel as 



K d (t,x) = 7 



1 



(47rf)5 



a il/2ad-x 



exp 



(x - 2-7T7I1) 2 

4r ' 



(A3) 



The expression in Eq. (|A1|) is fast convergent for large t, 
and the expression in Eq. (| A3|) is useful for small t. We 
can denote the truncated heat kernel K^(t,x) by 

K^(t,x) = Kd(t,x.) - Y exp(ir • x - tr 2 ). 

r£P d ,|r|<A 
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We apply the operator yi m {— «V X ) to heat kernels, 

£d,im(*> x ) = ^ m (-iV x )^(t,x). (A4) 

We can easily show that the zeta function has a rapidly 
convergent integral expression 



zLt^ 2 )= E 



3V(r) 
r 2 — g 2 



rGP d ,|r|<A 

+ (2,r)»jf * (e^ 2 /^ im (i,0)-^f) .(AS) 

To calculate the integrand, we use the Eq. (|A1|) when 
t > 1, and the Eq. (|A3|) in the case of i < 1. The cutoff 
A is chosen such that A 2 > Reg 2 . We can easily ver- 
ify that, when m\ = W2 (or equivalently a = 1), the 
Rummukaincn-Gottlicb's result in Ref. M is restored. 



Appendix B: The evaluation of the zeta function 

Z w {s;q 2 ) 

In this appendix we briefly discuss one useful method 
for numerical evaluation of zeta function -2io(s; q 2 )- Here 



we follow the methods and notations in Ref. |ll| . 

The definition of the zeta function Zf (s;q 2 ) in 
Eq. ([21 is 



Zf (s; q 2 )=J2 



r 3 



rEPd 



(r 2 — q 2 ) s 



(Bl) 



where the summation for r is carried out over the set 
P d = {rjr^r 1 (n+|d) , neZ 3 }, (B2) 

here the factor a is denoted in Eq. (|50"|) . The operation 
7 _1 is is defined in Eq. |QJ. We consider that the value 
q 2 can be a positive or negative. 

First we consider the case of q 2 > 0, and we can sepa- 
rate the summation in Zf (s;q 2 ) into two parts as 



E 

rGPd 



(r 2 — q 2 ) s 



E 



E 



r 3 



(r 2 - q 2 ) s ' ^ (r 2 - q 2 ) s ' 

(B3) 



The second term can be written in an integral form, 



E 



- y 



{ r 2 -q 2 Y T(s) , 

r z >q 2 r 2 >q 2 

rl 



At t s - l e- 1 ^-^ + I dt t^e-^ 2 -^ 



1 



r(«) Jo 



d^-VE^e-''"'- E 



'•:s 



(r 2 — q 2 ) s 



E ? ' 3 



(r 2 — q 2 ) s 



(B4) 



The second term neatly cancel out the first term in Eq. (|B3j) . Using Poisson's resummation formula we can rewrite 
the first term in Eq. (|B4[) as 



[\tt a - l e^ (-Yy n 3 sin(a7rn.d)e-^ n ) 2 /*, 



(B5) 



r 



where the imaginary parts are neatly canceled out. 

After gathering all terms we obtain the zeta function 
at s = 1 as, 

^ ZUl-q 2 ) = E^3 6 



V¥ y dte tq2 ^ 2 E n 3 sin(a7rn-d)e-^ n ) 2 / t . 



(B6) 



For the case of q 2 < 0, it is not necessary for us to 
separate the summation in Ziq(s; g 2 ), and it can be also 
written in an integral form. Following the same pro- 
cedures, we arrive at the same expression in Eq. (|B6I) . 



Hence, Eq. (|B6j) can be applied for both cases. 

Substituting d = (0,0, 1) into Eq. (|B6[) we obtain the 
zeta function in Eq. ([56| 



2 io(i;9 ) = 2^ 7,3 



r£P d 



r 2 — q 2 



\pK / die* 9 ^— J ri3 sin(a7rn3)e 



(B7) 



We can easily verify that, if mi = m<i, zeta function 
^io (1; <? 2 ) — > 0, the Rummukainen-Gottlieb's result is re- 
covered. 
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